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ABSTRACT 

Providing  precision  guidance  to  an  aircraft  landing  aboard 
a  ship  requires  a  robust  and  extremely  accurate  positioning 
system.  Several  carrier  controlled  approach  aids  (such  as 
the  AN/SPN-42A/46  Automatic  Carrier  Landing  System) 
currently  provide  the  required  accuracy  for  the  U.S.  Navy. 
However,  these  systems  experience  reduced  performance 
in  precipitation  and  are  difficult  to  maintain.  A  relative 
Global  Positioning  System  (GPS)  Kinematic  Carrier 
Phase  Tracking  (KCPT)  solution  has  the  potential  to 
provide  a  highly  accurate  solution  for  shipboard  landing 
operations  which  is  unaffected  by  weather.  There  are 
several  problems  with  a  GPS  based  shipboard  landing 
system  over  and  above  its  shore  based  counterpart.  The 
touchdown  point  and  the  GPS  reference  station  are  in 
motion  through  six  degrees  of  freedom,  the  ship’s 
dynamics  are  nearly  as  high  as  the  aircraft  dynamics 
(effectively  doubling  the  bandwidth  required  from  a 
navigation  sensor),  and  the  reference  station  experiences 
more  cycle  slips  and  masking  due  to  the  high 
electromagnetic  interference  environment  and  the  ship’s 
structure.  With  these  considerations  in  mind,  E-Systems, 
Montek  Division  has  developed  a  relative,  all-in-view, 
KCPT  solution  that  is  able  to  resolve  and  hold  cairier 
cycle  ambiguities  indefinitely.  Space  Vehicle  (SV) 
pseudorange  and  carrier  phase  data  as  well  as  ship’s 
motion  data  are  uplinked  to  the  landing  airaaft  to 
formulate  an  air  derived  solution.  The  system  features  a 
floating  point  solution  (ambiguities  not  resolved)  and 
employs  a  Kalman  filter  to  estimate  the  magnitude  and 
covariance  of  the  ambiguities.  The  results  of  the  Kalman 
filter  are  fed  to  an  ambiguity  resolution  algorithm  that 


operates  continuously.  The  output  of  the  GPS  solution 
(either  a  lloating  or  fixed  solution)  is  blended  with  high 
rate  inertial  data  in  a  filter  developed  by  the  Naval  Air 
Waifaie  Center  Aircraft  Division  (NAWCAD),  Patuxent 
River  to  provide  a  final  position  solution.  This  three  step 
process  provides  a  continuous  navigation  solution  through 
loss  and  reacquisition  of  SVs  and  cycle  slips.  The  system 
has  been  evaluated  with  two  sets  of  approach  and  landing 
data.  Tlie  first  set  was  collected  during  an  FAA  Category 
III  Feasibility  Study  in  which  E-Systems  successfully 
completed  100  approaches  and  landings  with  a  West  wind 
1124  aiiplane  using  laser  tracker  to  provide  ‘truth’ 
position  data.  The  second  set  was  collected  during  a  series 
of  at-sea  approaches  by  a  GPS  equipped,  Naval  Rotary 
Wing  Test  Squadron,  SH-60F  Seahawk  helicopter  to  the 
aircralf  carrier  USS  ENTERPRISE  (CVN-65).  The 
production  GPS  antenna  (used  for  these  tests)  on  the  SH- 
60F  is  located  in  an  area  susceptible  to  SV  masking  and 
cycle  slips,  and  the  data  collected  during  the  approach 
profiles  contained  multiple  occurances  of  both.  In  all  cases 
the  KCPT  solution  provided  a  seamless  fixed  ambiguity 
approach  solution  with  an  average  of  42  seconds 
ambiguity  fixing  time. 

INTRODUCTION 

Probably  the  most  challenging  task  in  aviation  is  landing 
abiwd  mi  aircrrft  carrier.  Night  operations,  foul  weatha, 
and  pitching  decks  do  nothing  to  help  the  situation. 
Howevo'  the  ability  to  laundi  and  recover  aircraft  in  these 
conditions  adds  to  the  carrier's  potency.  The  U.S.  Navy 
cumenlly  uses  several  systems  to  aid  the  pilot  during 
landing.  The  AN/SPN46  is  the  most  capable  of  these. 
Using  a  shipboard  precision  approach  radar  and  racfio  data 
link,  the  system  provicbs  both  guidance  needles  and 
automiiic  control  to  aircraft  equipped  with  a  radar  beacon 
transponder  and  data  link  receive.  The  AN/SPN41  is  an 
pulse-aule  scanning  beam  system  that  provicbs  civilian 
Instmnent  Landing  System  (ILS)  look  alike  needles  and  is 
used  as  m  independent  monitor  to  the  AN/SPN46.  In-close 
optical  aids  and  a  Landing  Signal  Officer  (LSO)  who 
monitoi's  the  approach,  complete  the  suite  of  aids,  and  are 
used  during  all  approaches.  Althoi^h  having  supported 
thousaids  of  successftil  lancfings,  the  AN/SPN46  suffeis 
some  shoita^mings.  Only  two  aircraft  can  be  provid'd 
guichncc  simultaneously,  a  transponder  is  required  in  the 
aircrrft  to  iK'hieve  the  required  accuracy,  performance  can  be 
Irmital  in  heavy  rain,  and  the  system  is  relatively  complex 
to  maintan. 

Currently  the  naval  aviator  must  rely  on  Precision 
Appnxch  Radar  (PAR)  or  nonfrrecision  approa:h  aids  such 
as  TAG  AN  at  most  civilian  and  U.S.  Air  Force 
installations.  With  civilian  interest  in  GPS  for  shore  based 
approa'hes  and  advances  in  GPS  processing  techniques,  the 
Navy  saw  opportunities  to  increase  shore  based  bad  weather 


landing  capability,  reduce  the  amount  of  ecjuipment  required 
for  lancing  operations,  and  provicb  a  consistent  set  of 
procediies  andcispl^s  for  both  shore  based  and  shipboard 
operations.  Future  compatibility  with  civil  GPS  landing 
systems  is  also  possiWe.  However,  the  shipboard 
environment  presorts  unique  challmges  ova*  and  above 
those  of  landing  ashore  The  aircraft  must  land  in  a  zone 
approximately  20  meters  long  and  3  meters  wide  (as 
opposed  to  the  several  hundred  meter  lancing  box  shore 
based^  which  requires  decimeta  accuracy  from  the 
navigation  sensor.  In  additbn,  the  ship  (and  thus  the  GPS 
refaoice  station)  translates  and  rotates  in  six  degrees  of 
freedom.  Multipath  effects  on  the  refaaace  statiai  and  SV 
block^e  is  of  special  concern  since  the  environment  aboard 
ship  is  (fynamic,  and  siting  relative  to  otha  stationary 
objects  cannot  be  contrciled  The  AN/SPN46  uses  well 
developed  tedmiques  to  compensate  for  the  effects  of  ship’s 
motion  and  to  provicb  control  and  dspl^s  to  the  aircraft 
[1],  so  the  first  aspect  addressed  was  providng  a  robust, 
accurate,  relative  position  solution  at  a  high  enough  sample 
rate  to  feed  these  stabilization  and  control  algorithms. 

A  carrier  phase  ambiguity  resol\ed  GPS  solutbn  is  the 
basis  of  the  system  shown  herein  Differenced  acceleration 
measurements  from  the  aircnft  and  ship  mitigate  the  low 
(1-2  Hz)  data  rate  typicri  of  GPS.  GPS  pseudorange  and 
carrier  phase  data  from  the  reference  station  as  well  as  ship’s 
inertial  data  are  up  linked  to  the  lan(ing  aircraft  to  formulate 
an  air  derived  relative  position  solution  at  10  Hz.  The  ship's 
GPS  (fata  is  up  linked  at  1  Hz  and  the  ship’s  inertial  data  at 
10  Hz.  The  GPS  processing  algorithm  provicbs  a 
continuous  floating  point  ambiguity  solutbn.  This 
solution  has  accuracies  typical  of  carrio*  smoothed  codb 
differential  GPS  (1-2  metas)  andprovi(bs  adequate  accuracy 
for  aircraft  contrd  at  range.  Once  ambigiities  are  resolved, 
a  fixed  solution  of  higher  accuracy  is  provi(bd  This 
solution  has  the  recjuired  accuracy  for  aircraft  control  to 
touchcbwn  (less  than  30  cm).  Satellite  switching  during  the 
ambiguity  search  process  is  addressed  so  outages  are  not 
experiaiced  as  long  as  5  common  satellites  are  maintained 
in  view.  The  relative  GPS  solution  and  the  d&ffeienced 
acceleration  measurements  are  used  in  a  blard  filta  to 
provicb  a  high  rate  positbn,  velocity,  and  acceleration 
solution  to  the  stabilization  and  contrd  equatbns.  This 
blend  filter  uses  a  Kalman  filta  to  estimate  inertial 
measurement  biases  and  to  mitigate  the  latency  of  the  GPS 
solution  (typically  one  to  two  secondi).  With  better 
estimaes  of  the  inertial  measurements,  the  navigation 
solution  can  be  coasted  through  period;  of  GPS  outages, 
adding  to  continuity  of  service. 

The  design  of  the  KCPT  software  and  inertial  blend  filter  is 
presented  along  with  results  from  two  sets  of  data  The  first 
set  is  several  ^proaches  of  a  Westwind  1124  turbojet  to  a 


fixed  base.  The  second  set  is  shipboard  approaches  of  an 
SH-60F  heliccpter  to  the  USS  ENTERPRISE  (CVN-65). 

KINEMATIC  CARRIER  PHASE  TRACKING. 

This  paper  describes  a  technology  that  calculates  a  highly 
accunte  positbn  solution  that  can  survive  in  an  aircraft 
carrier  environment.  The  main  component  of  this 
technology  is  the  GPS.  GPS  is  a  satellite  based navigftion 
system  that  continuously  transmits  timing  frequaicy,  and 
SV  position  information  on  the  LI  and  L2  channds  to 
potential  users.  The  full  constdlation  consists  of  24  S  Vs  in 
half  geosyichronous  orbits.  The  constdlation  is  contrdled 
by  the  Department  of  Defense  (DoD)  which  monitors  the 
position  and  clock  accuracy  of  the  satellites.  The  timing 
information  is  a  coded  signal  that  allows  the  user  to 
detamine  the  time  elapsed  for  the  signal  to  transvCTse  the 
(istance  between  the  S  V  and  the  user.  By  knowing  the  time 
the  signal  left  the  SV,  and  the  speedof  the  signal  (speed  of 
lighti  the  user  can  determine  the  distance  from  itself  to  the 
SV.  This  range  measurement  is  biased  by  clock 
inaccuracies  in  the  SV  and  user  clocks  Because  of  these 
timing  errors  it  is  not  a  true  range,  but  a  “psuedrrange.”  By 
knowing  the  position  of  the  SV  (ephemeris  data),  and  die 
distance  from  itself  to  the  S  V,  the  user  can  triangulate  its 
own  position. 

The  constant  motion  of  the  SVs  and  possible  user  motion 
dictate  that  the  receiver  must  also  be  able  to  track  the 
change  in  frecjuQicy,  or  Doppler  shift.  Integration  of  die 
Dopplff  shift  over  time  yields  a  highly  accurae  delta  range 
measurement  which  is  proportional  to  the  advance  in 
signal-carrier  phase  betweai  two  time  epochs.  Taking 
advantage  of  these  techniques  typicdly  provirfes  a  delta 
range  measurement  to  an  accuracy  of  a  tenth  of  the 
wavelength.  LI  wavelength  is  19  cm  while  L2  wavelength 
is  24  cm  which  yields  an  deltamnge  accuracy  of  1-3  cm. 
Tmc  range  between  the  SV  and  receiver  cannot  be 
determined  by  integrated  Doppler  techniques  alone  becau.se 
the  eonstait  of  integration  is  ambiguous.  Technicjres  will 
be  (iscu.ssed  on  how  these  whole  number  ambiguities  can 
be  resolved  The  interferometric  concept  of  using  integrated 
Doppler  to  resolve  a  relative  veetor  betweai  two  antennas 
will  be  referred  to  as  KCPT  and  involves  five  steps:  1) 
Tracking  the  carrier  phase  advance  of  the  SV  carria-  signal, 
2)  Perfonning  a  double  difterence  on  the  raw  carrier  phase 
and  pscudnrange  measurements,  3)  Solving  for  the  relative 
vector  between  the  ground  antenna  and  airborne  antenna, 
using  double  differenced  pseudnrange  measurements,  4) 
Calcukling  a  highly  accurate  relative  vector  between  the  air 
and  ground  antennas  by  resolving  carrier  phase  ambiguiiie.s, 
and  5)  Compensating  for  ship’s  motion  and  high  frecfiaicy 
airatft-ship  kinemaics  with  ship  and  air  inertial  data. 

Carrier  Phase  Tracking.  Once  a  receiver  channel  locks  on 
to  a  cairier  signal  from  a  single  SV,  the  channd  keeps  a 


mnning  count  of  the  cycles  based  on  the  Doppler  shift 
present  on  LI  and  L2.  This  is  done  by  integrating  the 
Doppler  shift  ovct  the  interval  of  the  qioch.  At  the 
(inclusion  of  eadi  qioch,  the  estimae  of  the  caniff  phase 
count  of  that  epocii  is  added  to  the  count  of  the  previois 
epoch  to  keep  a  runnir^  total  of  the  number  of  carriff  phase 
cycles.  As  long  as  the  reedva"  keqis  lode  on  the  carria,  the 
deltarange  measuremert  is  extremely  accxrrae.  The  receiver 
must  begin  to  “counf’  or  integrate  at  some  point.  When  the 
receivo'  begins  this  count  or  if  the  receiver  looses  lock  on 
the  carria,  it  has  no  knowledge  of  the  previars  Dopplff 
shift  count,  so  it  begins  to  integrate  at  some  arbitrary 
whole  number.  Therefore,  the  true  whole  number  is 
ambigtxrus,  and  must  be  resolved  to  obtain  an  accurae 
position  fix  using  carrier  phase  measuronents.  The 
methock  of  fixing  ambiguities  will  be  discussed  later. 

The  convetiional  goal  for  the  KCPT  process  is  to  be  able 
to  track  LI  frequencies  and  resolve  LI  ambigirities. 
However,  LI  ambiguities  are  difticult  to  resolve  because  of 
the  small  LI  wavelaigth.  Differencing  LI  and  L2  earner 
phase  measurements  will  produce  a  wide  lane  wavelaigth  of 
86  cm.  The  penalty  for  using  the  wide  lane  measurement  is 
that  the  noise  level  is  inaeased  over  the  LI  measurement 
by  a  factor  of  six.  Certainly,  the  high  accuracy  is  obtained 
when  LI  carrier  cycle  is  used,  but  the  wide  lane 
measurement  can  prove  as  a  very  effective  intermeeSate  step 
in  resolving  LI  ambiguities.  The  Ashtedi  Z-12  receiver  has 
proven  effective  in  tracking  LI  andL2  without  cartia  phase 
cycle  slips  that  adversely  effect  operations.  The  receiver 
employs  a  Z-Code  tracking  capability  that  ptovi±s  high 
quality  code  and  carrier  phase  measurements  on  both  LI  and 
L2  frequaides  even  with  encrypted  timing  codes. 


Double  Diffaence  Processing.  Perfonning  a  double 
difference  on  the  raw  amier  phase  measuements 
accomplishes  two  important  tasks.  It  processes  the  range 
measurements  so  they  can  be  rearfily  applied  to  a  relative 
vector  solution,  and  cancels  out  most  of  the  systematic 
errors.  The  KCPT  solution  solves  for  a  relative  vector 
between  the  ground  and  airborne  antennas.  Therefore,  the 


range  measuiements  of  the  same  S  V  must  be  diffeiEnced 
betwem  the  two  receivers  to  obtain  the  difference  in  range 
measurements.  Referring  to  Figure  1,  two  range 
measuiements  are  ciffeienced  for  SVi.  The  vector  betweoi 
the  two  receivers  is  called  lb.  By  differencing  the  range 
measurement,  the  result  is  the  scaler  bi,  which  is  the  dot 
product  of  the  baseline  vector  and  the  unit  direction  vector 
(ui).  The  unit  direction  veaois  indicate  the  diecdon  of  the 
SVs  relative  to  a  point  halfw£^  betwem  the  ship  and 
airborne  positbn.  This  single  dfiference  is  performed  for  all 
SVs  in  view.  Ui  is  assumed  equal  for  the  two  receivers 
because  the  dfference  in  unit  vectors  from  the  ground  to 
S  V  and  airborne  to  SV  is  negligible.  This  single  dfference 
virtually  cancels  out  common  errors  between  ground  and  air 
causedby  the  SV,  Selective  Availi)ility  (SA),  and  most  of 
the  atmospheric  mors. 

The  double  dffamce  is  performed  by  dfifoencing  the  single 
dfference  of  a  target  SV  against  all  other  SVs  in  view.  This 
target  SV  is  usually  the  highest  elevation  SV  for  severd 
reasons.  This  choice  usually  yidds  the  best  geometry  and 
tropojpheric  error  mitigation,  and  the  highest  elevation  SV 
will  remain  in  view  the  longest,  thus  assuring  the  target 
SV  will  not  have  to  be  switched  for  visibility  reasons 
during  an  approadi.  The  advantage  of  performing  the  double 
dfference  is  that  the  receiver  clodc  mors  cancel  Referring 
to  Figure  1,  it  can  be  sem  that  the  double  dffamce  is  the 
dot  produa  of  tme  relative  vector  with  the  dfference  of  the 
two  unit  drection  vectors,  hi  the  case  of  pseudoranges,  the 
double  diffeience  yields  a  relative  range  measurement  with 
no  ambiguities  but  significantly  larger  noise  than  the  phase 
measuiements.  The  carrier  phase  double  dfference 
measuiements  yields  a  low  noise  relative  range 
measuiement  but  with  double  dflfe^nce  ambiguities. 

Pseudrange  Double  Diffeience  Solution.  The  goal  for  this 
process  is  to  produce  a  relative  vector  solution  based  on 
canier  phase  measuiemmts.  Since  the  carrier  phase 
ambiguities  are  still  unknown,  a  relative  position  solution 
is  not  possiWe  until  the  ambiguities  are  calculated  An 
initial  estimate  of  the  integer  ambiguities  can  be  performed 
by  comparing  the  double  dffeienced  carrier  phase 
measuiements  with  the  double  dffeienced  pseudorange 
measuiements.  Referring  to  Figure  1,  an  estimate  of  this 
integer  ambiguity  can  be  made  by  Equation  1.  Since  typical 
smoothed  pseudorange  double  dfference  error  has  a  3a 
value  of  2  m,  and  an  LI  wavelength  is  19  cm,  the  LI 
integer  ambiguity  can  be  in  error!  11  wavelengths. 

ddpri2  ~  ddcpn 

N, - (1) 

where 

Ni  =  initid  estimate  of  double  dfifaence  integer 
^unbiguity. 

=  pseudorange  double  dfference  betweai  S  V,  and 


SV2. 

dd(,p  =  carrier  phase  double  dfiference  between  SVI 
andSV2. 

X  =  waveloigth. 

The  Pseudorange  double  dfference  can  also  be  used  to 
calculate  a  relative  vector.  This  is  advantageous  because  it 
can  serve  as  a  backup  to  the  carrier  phase  solutm  if  carrier 
phase  measuTonents  are  not  available.  The  double  dfference 
code  solution  can  also  configure  the  algorithms  that  evolve 
into  a  KCPT  solution.  By  solving  for  the  b  vector  in 
Eolation  2,  an  estimate  of  the  relative  veaor  can  be 
produced  which  should  meet  the  shipboard  accuracy 
requirements  for  a  manual  approxh  to  1/2  nmi. 

Resolve  Carrier  Phase  Ambiguities.  The  LI  carrier  phase 
measujEments  are  highly  accurie  measurements  that  can 
give  a  relative  position  solution  to  centimeter  level 
accuracy.  To  accomplish  this,  the  double  dfference 
ambiguities  must  be  resolved  Equatbn  2  is  expanchd  in 
ecjiation  3  to  inclucfe  double  dfference  ambiguities.  Note 
that  this  process  is  first  initialized  with  the  ambiguities 
first  estimated  in  equation  1.  Once  the  ambiguities  have 
been  initialized  they  can  be  in  error  (typically  less  than  ±11 
wavelengths).  This  yields  23“  possiWe  combinations.  There 
are  several  method  for  resolving  these  ambiguities,  three 
categories  will  be  dscussed  1)  Numerical  search 
algorithms  2)  Ground  based  carrier  phase  measuiements  that 
force  a  high  rate  of  ciiange  of  geometry,  and  3)  a  Kalman 
filter  that  estimates  the  ambiguities  and  covariance  of 
ambiguities. 
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where 


Njj  =  initial  estimie  of  cbuble  dfferaice  intega 
ambiguity  between  SVI  and  S  V  j. 
d4pij=  carrier  phase  cbuble  dfference  between  S  v^  and 
SV, 

X  =  initial  ambiguity  estimate. 

Ujj  =  dfference  between  unit  veaors  of  S  V  i  andj. 
n  =  number  of  SVs  used  in  the  solution. 


The  first  method  to  iesol\e  carrier  phase  ambiguities  is  a 
numerical  method  using  a  algorithm  search  methocblogy 
that  tries  dffeient  combinations  of  ambiguities  and 
examines  the  residuals  [2].  Residuals  near  zm)  are  identified 
and  kept  as  a  candichte  ambigdty  sets.  The  variance  of 
these  residuals  arc  be  monitorcd  to  determine  the  correct 
ambiguity  set. 

Another  method  of  resolving  ambiguities  is  to  provitb  a 
ground  based  carriff  phase  measuiement  bubble  that  the 
aircraft  will  be  required  to  fly  ova  [3].  With  the  resulting 
fast  changing  geometry,  the  ambiguities  become  readly 
observable,  and  the  extensive  search  algorithmspiesented  in 
the  previdis  paragiaph  arc  not  required 

The  third  methodemploys  a  Kalman  Filter  to  estimate  the 
amhiguiiies.  The  K^man  filter  uses  differenrcs  in 
psuedmange  and  carrier  phase  double  difference 
measurements  as  well  as  aircraft  (tynamics  to  converge  on 
the  ambiguities.  A  side  benefit  of  this  method  is  ajow 
noise  md  accurate  “intennediate  solution”  known  as  the 
float  solution,  that  can  be  used  while  the  systan  is 
resolving  ambigiities.  The  Kalman  filta  alone  is  not  as 
reliable  as  the  numerical  estimstion  approach  because  the 
Kalman  only  calculates  the  diffeaence  betweai  the 
psuedtTiuige  and  carrier  phase  measurements.  Other  errors 
caused  by  the  atmosphae  and  geometry  arc  seen  by  the 
Kalmai  filter  as  part  of  the  ambiguity.  Therefore,  the 
Kalimui  filler  must  be  augmented  by  a  method  that 
optimizes  the  three  dmensional  geometry  of  a  position 
solution  into  the  geometry  of  the  current  SV  constellation. 
An  effective  geometry  optimization  method  is  known  as 
the  Teunissen  method [4]. 

Inert kil  Blending.  To  provid  output  betweoi  the  low  rate 
(typically  1  to  2  Hz)  GPS  measurements,  acceleration 
measuiements  from  the  aircraft  and  ship  are  used  to 
extnipolatc  from  the  last  valid  position  determined  from 
GPS  [5].  Using  this  technique,  the  stabilization  and  control 
laws  am  be  given  10  to  20  Hz  data  whidi  is  adequate  to 
tmek  the  cx)mbined  ship  and  aircraft  motioa  Two 
conseaitive  GPS  position  solutions  are  used  to  determine 
an  initiri  velocity,  and  the  second  GPS  solution  is  used  as 
an  initiii  position.  From  these,  the  acceleration 
measuiements  arc  cfifferenced  and  integiated  (assuming 
constiiii  acceleration  ova  the  sample  time)  for  velocity  and 
position  pmpagation. 

Tlie  transformation  from  the  Earth-Caitaed-Earth-Fixed 
(ECEF)  frame  to  the  East-Nbrth-Up  (ENU)  frame  is 
includxl  in  the  GPS  processing,  so  the  blend  filta  is 
presented  with  an  ENU  referaiced  relative  position  vector 
from  the  ship’s  GPS  antenna  to  the  aircraft’s  GPS  antenrtL 
If  dilTeient  refaence  points  are  desired  (for  example  from  the 
ship's  touchdown  area  to  the  airosff  s  tailhook),  moment 


arm  corrections  for  these  offsets  arc  done  in  the  ENU  frame. 
Thus,  the  offsets  whidi  are  known  in  the  body  frame  must 
be  transformed  to  the  ENU  frame.  This  transformation  uses 
the  platform’s  Euler  attitude  measurements  and  is  simply 
the  produa  of  three  rotations;  roll  first,  then  pitch,  then 
yaw  [5]. 

p  =  rx  (4) 

Where: 

r  =  is  the  standard  Eula  angle  conversion  matrix. 

P  =  is  the  ENU  oriented  vector, 
x  =  in  the  body  axis  vector 

Using  this  rotation  matrix,  two  moment  arm  correaions 
arc  applied  to  the  relative  GPS  solution;  the  moment  arm 
from  the  ship's  GPS  antenna  to  Inertial  Measurement  Unit 
(IMU)  and  from  the  aircraft’s  GPS  antoma  to  MU.  This 
allows  comparison  of  the  relative  GPS  solution  with 
acceleration  measuiements  in  subseqient  processing 
without  referaidng  the  acceleration  measurements  to  otha 
points  on  the  platfonn  (whidi  requires  knowlajge  of 
attitude  rate  and  attitude  acceloation).  If  a  stiapdbwn  sensor 
is  used  for  measuring  acceleration,  the  acceleration 
measuiements  arc  transformed  to  the  ENU  frame  in  the 
propagation  of  the  GPS  solution  using  the  same  rotation 
matrix  above.  If  the  accelerations  arc  taken  with  a  gimbaled 
sensor;  this  transformation  would  not  be  required  The 
following  equations  arc  used  to  propagate  betweai  GPS 
updates  with  body  axis  acceleration  measuiements. 
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Where 

dt  =  time  since  the  last  propagation. 

F  =  body  to  ENU  transformation  matrix. 

ac  =  aircraft, 

s  =  ship. 

P  =  East  North  Up  relative  position  vector. 

A  =  body  axis  acceleration  vector, 

k  =  current  epoch. 

k-1  =  previous  epoch. 

I  =  3x3  Identity  matrix. 

0  =  3x3  zero  matrix. 

Given  two  perfed  relative  GPS  solutions  to  initialize  the 
propagation  and  perfed  acceleration  measuiements,  no 
further  GPS  measuronents  would  be  required  Since  this  is 
obviously  not  the  case,  some  means  of  continuously 
updating  the  inertial  propagation  is  needed  This  process 
should  dstir^uish  between  enors  in  the  initial  conefitions 
for  position  and  velocity  and  acceleration  measuiement 
errors.  A  Kalman  filta  is  used  to  compaie  the  position 
output  of  the  inertial  propagition  with  the  relative  GPS 


solution  to  make  this  allocation.  As  an  added  benefit,  the 
Kalman  filto-  inclucfes  the  effects  of  GPS  measuiement 
latency,  so  this  latency  is  mitigated  The  outputs  of  the 
Kalm^ui  filter  are  correaions  to  the  cunent  position  and 
velocity  in  the  inertial  propagation  (acjusting  the  initial 
condiiions  of  future  integiations)  and  corrections  to  all 
future  ijcceleration  measuiements  (the  acceleration 
measuiement  bias).  The  only  buffering  requiiement  is  that 
the  old  positbn  outputs  from  the  inertial  propagation  be 
saved  so  the  output  matching  the  GPS  measuiement’s  time 
can  be  used  in  forming  the  Kalman  measuiement  vector. 

The  Kidman  filter  models  position  and  velocity  initial 
condition  errors  as  constant.  Body  axis  acceleration 
me<asuiement  errors  are  also  modeled  as  constait.  These  are 
not  bad  assumptions  since  the  initial  condition  errors  are 
constaii,  iuid  for  short  run  times  driring  an  approach,  the 
acceleration  measuiement  errors  are  nearly  constant. 
Howevu*,  for  this  filter  the  Kalm^  statistics  for 
accelantion  bias  are  modeled  as  integmted  white  nois^  (or 
random  walk). 


The  intent  here  is  not  to  rederive  the  Kalman  filter 
equati(^ns,  so  the  measuiement  and  state  equations, 
(equmions  6-10)  are  given  in  the  standard  form  without 
furtha  exphuiation. 


Zk  =  Hk  Xk  +  Vk 

Xk+j=  ^kXk-^GkWk 
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Pj^  =  initid  condtion  position  error. 

Ahi,,s=  acceleration  bias  vector  in  body  axis. 

dT  =  time  since  last  Kalman  step  plus  GPS  latency. 

G  =  noise  transfer  function. 

<I>  =  zero  if  corrections  are  fed  back. 


The  measuiement  vector  (AE,  AN,  AU)  is  the  diffeience  in 
the  ENU  output  from  the  inertial  propagation  and  the 
relative  GPS  solutbn.  Thus  the  filter  is  opaaling  on  enror 
terms  and  the  outputs  are  inaemaital  correaions.  These 
correaions  can  either  be  fed  back  to  the  inertial  propagation 
(in  which  case  the  Kalman  output  should  go  to  za-os)  or 
can  be  accumulated  and  ^plied  to  the  output  of  the  inertial 
propagation  (in  which  case  acceleration  biases  cause  linear 
growth  in  velocity  correaions  and  exponaatial  growth  in 
position  correaions).  To  reduce  the  chance  of  ovaflowing 
computer  representations  of  these  correaions,  the  outputs 
for  positbn  and  velocity  correaions  are  fed  back  to  the 
inertial  propagation,  and  outputs  for  acceleration  correaions 
are  accumulated  and  applied  to  subseqient  raw  acceleration 
measuiements. 

Terms  in  the  measuiement  noise  covariance  matrix  and  the 
process  noise  covariance  matrix  are  determined  by  the  noise 
on  the  relative  GPS  solutbn  and  on  how  bias-like  are  the 
errors  on  acceleration  measuiements.  Since  the  relative 
GPS  solutbn  is  expected  to  have  very  low  noise  (any  error 
is  likely  to  be  bias-like  which  the  filter  does  not  modeU 
the  diagonal  terms  of  the  R  matrix  (E[vv^)  are  relatively 
small  (typically  5  cm  or  less).  The  nonzero  cfiagonal  toms 
of  the  Q  matrix  (Efww”^])  have  been  chosen  empirically  at 
betweai  0.2  and  0.5  mVs"^  to  provicfe  reasonable  estimrtion 
of  acceleration  measuiement  biases  while  still  tracking 
dynamics  (all  non-diagonal  tenns  in  the  Q  matrix  are  set  to 
zero).  The  relatbnship  betweai  the  acceleration  bias  state 
process  noise  and  the  position  and  velcKity  initid  condtion 
noise  is  in  the  G  matrix.  So  the  terms  in  the  Q  matrix 
corresponding  to  the  positbn  and  velocity  error  states  are 
zero.  The  G  matrix  simply  inducts  the  effects  of 
integiating  the  random  walk  process  once  for  velodly  noise 
and  again  for  positbn  noise  [5]. 

Figure  2  shows  the  relationship  betweai  the  inertial 
propagation  and  the  Kalman  filter  in  the  feedback 
configuration.  PVA  denotes  position,  velocity,  and 
acceleration  respeaively.  The  'u’  subscript  means 
unconected,  and  the  'c'  subscript  means  correaed  Thus  Pc, 
Vc,  and  Ac  are  the  outputs  that  are  used  in  the  stabilization 
and  contrci  algorithms. 
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Figure  2  Inertial  Propagation  and  Kalman  Blend  Filter 


The  Kalman  gain  and  state  covariance  matrices  are 
pnopagited  in  time  with  the  standard  equations  [5]  [6].  The 
outputs  of  the  inertial  propagation  should  trade  the  vehicle 
(fynamics  during  the  approach,  and  givai  small  latency  on 
the  acceleration  measurements  will  providb  currerit  relative 
position,  velocity,  and  acceleration  between  the  platfam 
IMU’s  in  the  ENU  frame  to  the  stabilization  and  contrd 
equations.  To  correa  for  offsets  betweaa  the  IMU's  and  the 
desiroJ  icicTOice  points,  moment  arm  correaions  can  be 
applied  to  the  output  of  the  blendfilter.  The  position  vector 
reference  endpoints  can  be  moved  using  only  the  moment 
aim  veciois  and  the  vehicle  attitudes.  Howevo-,  to  move  the 
velocity  neferoice  requires  attitude  rate  information,  and  to 
move  the  acceleration  requires  attitude  rate  and  attitude 
accelantion  information,  which  up  to  this  point  in  tlie 
processing  are  not  required. 

Once  the  outputs  of  the  blend  filter  are  reterenced  to  the 
desiiui  points  on  the  ship  and  aircraft,  similar  stabiiizcttion 
and  control  equatbns  as  used  in  the  AN/SPN46  are  applitri 
[1].  This  includbs  ship’s  yaw  notch  filtering  and  low  pass 
filtering  before  use  in  the  ENU  to  body  rotatbn.  Ship's 
heave  compensation  will  also  be  calculated  so  the  aircraft's 
position  is  not  oontrdled  to  match  the  ship’s  vertical 
motion  until  the  last  12  second  or  so  of  the  approach. 


methodblogy,  triple  redunchney,  and  indepaidence  of  the 
final  solution  from  the  ambiguity  fixing  process.  The 
design  also  allows  integrity  checking  at  several  points. 

Observation  Matrix  Calculation.  For  the  RKCPT  solution, 
the  positbn  of  the  base  antenna  is  not  known,  therefore  the 
unit  direction  vectors  are  calculated  by  an  ita:ative  least 
sejuares  method  The  ground  unit  cfirection  vectors  are 
calculated  and  averaged  with  the  airborne  unit  direction 
vectors  to  obtain  the  solutbn  unit  drection  vectors. 

Kalman  Filter  Implementation.  The  Kalman  filter  is  the 
primary  tool  in  estimating  the  cania*  phase  ambiguities. 
The  Kalman  filter  estimates  the  aircraft  dynamics  as  well  as 
the  floating  point  wide  lane  ambigiities.  Note  that  this 
implementation  estimates  the  wide  lane  ambiguities.  This 
is  done  for  two  primary  reasons.  The  first  being  that  the 
reejuired  ebtn  up  linked  to  the  aircraft  from  the  ship  is 
reduced  Also,  wide  lane  ambiguities  are  more  reacily 
calculated  than  LI  ambiguities,  and  wide  lane  accuracies 
should  be  more  than  adequate  for  shipboard  approach  and 
landng. 

The  Kalman  filter  has  8  +  n  states:  Three  states  each  for 
position,  velocity,  and  acceleration  plus  n-1  double 
difference  ambiguities.  The  v^uiable  n  is  the  number  of  SVs 
in  view.  This  is  a  straight  forward  implementation  except 
for  the  fact  that  the  constdlaiion  is  not  static  The  Kalman 
filter  must  be  able  to  ch^uige  size  from  epoch  to  epoch 
because  of  SV  masking  and  cycle  slippage.  Therefore,  each 
ambiguity  state  must  be  evaluaedeach  second  to  detamine 
one  of  three  outcomes:  1)  The  ambiguity  state  will  retain 
the  covariance  information  as  well  as  the  current  floating 
point  estimate  of  the  ambiguity,  2)  The  ambiguity  state 
will  reset  because  of  a  cycle  slip  or  the  process  has 
detennined  that  the  previous  initialization  of  the  ambiguity 
was  insufficient,  and  3)  The  Kalman  filter  will  drop  or  add 
the  ambigiity  to  it’s  states.  The  accorxion  nature  of  this 
Kalman  filter  is  a  key  to  the  robust  performance  of  the 
RKCPT  process. 


RELATIVE  KCPT  DESIGN. 

Tlie  relative  KCPT  (RKCPT)  employs  the  same  basic 
concepis  as  the  standard  KCPT  solution.  Howevo-,  there  are 
signiliciuil  differences  that  require  an  independent  design  of 
the  KCPT  solution.  The  first  difference  is  that  the 
touch(t)wn  point  is  translating  in  three  degrees  of  freedom 
and  the  nmway  is  rotating  in  three  degrees  of  freedom.  The 
second  dlfcience  is  that  the  aircrrft  carrier  envircTunent  is 
hostile  to  the  KCPT  process.  Multijie  SV  maskii^s  md 
canid-  phase  cycle  slips  require  a  robust  method  to 
accommodtie  the  relative  fast  changing  SV  constdlaiion. 
The  five  design  features  that  are  unique  from  the  standrd 
KCPT  solution  are  calculating  the  observation  matrix, 
implementation  of  the  Kalman  Filter,  ambiguity  fixing 


Ambigiity  Fixing  Methodrlogv.  The  RKCPT  method  uses 
the  Teunissen  method  to  resolve  the  ambiguities.  The 
ambiguity  estimaor  desaibed  by  P.  J.  G.  Teunissen 
consists  of  applying  a  matrix  transform  (known  as  the 
Gauss  Transform)  to  a  set  of  double  differenced  floating 
point  ambigiities  with  the  aim  of  dea-)upling  their 
statistical  uncaiainties.  If  the  algorithm  successfully 
reduces  the  covariances  -  not  the  pure  variances  -  of  the 
ambiguities,  the  “integer  least  squares”  estimetion  of  the 
integer  ambiguities  becomes  much  simpler  and  presents  the 
system  with  statistiailly  reliable  integer  wavelorgth  counts. 

For  simplicity  consitbr  a  system  of  just  two  integer 
ambiguities.  If  a,  and  a2  represent  the  original  estimite  of 


the  double  diffaence  ambiguities,  then  the  application  of 
the  simple  Gauss  transfonn, 

r  1  01 

1  01) 
a  1 

produces  the  ambiguities  and  ^  with  explicit 
reprcsentaticm  given  by 

(12) 

(13) 

The  algorithm  proposed  by  Teunissen  ^plies  “integer  least 
squmes”  estimation  to  this  linear  combinations  of  the 
originid  tunbiguities.  This  operation  proves  successful  if 
the  inuislonned  ambiguities  ai  anda2  possess  less  statistical 
uncertainty,  as  indicated  by  their  covariances,  than  the 
origin^d  iunbigiities,  aj  anda2.  In  that  case,  the  problem  of 
estimating  the  transformed  ambiguities  is  more 
straightforwaid  The  Gauss  transfonn  uses  the  more  cmain 
ambiguities  to  assist  in  the  determination  of  the  less  cotain 
ambiguities.  Differing  geometries  of  the  S  Vs  relative  to  the 
base  receiver  prodhee  more  or  less  uncertainty  in  the 
estijnaion  of  the  double  diffeience  ambiguities.  To  retrieve 
the  cstiinjtes  of  aj  and  a2  after  obtained  the  integer  estimates 
of  a,  tUid  we  simply  apply  the  inverse  of  the  Gauss 
transfonn. 

The  whole  process  compiises  a  kind  of  “geneicdized’ 
rounding  algorithm  producing  integer  ambiguity  estimaes 
from  the  input  floating  point  ambiguities  in  such  a  way  as 
to  reduce  their  uncertainties.  In  many  cases  the  algtrthm 
may  output  some  of  the  same  integer  ambiguities  as  would 
a  simple  m unding  operation  perfoimed  on  the  origirtd 
lloaling  point  ambiguities,  but  even  in  those  cases,  we 
wouldhave  inaieased statistiad  amficbnce  in  the  results  by 
perldiming  the  Teunissen  procedire.  Figure  3  desciibes  in 
general  lerins  the  Teunissen  method 
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Figure  3  Generalized  Teunissen  Method 


Triple  Redunchnev.  The  KCPT  process  requiies  a  taiget  S  V 
to  ciffeience  all  otho*  SVs  against  in  the  double  diffeience. 
If  a  cycle  slip  or  masking  occurs  on  any  other  S  V  than  the 
target  SV,  the  process  will  simply  reset  that  particular 
[imbiguily,  and  the  KCPT  fixed  solution  c^in  simply 


continue  without  the  S  V  in  question.  But  if  a  cycle  slip  or 
masking  would  to  occur  on  the  target  SV,  then  the  whole 
ambiguity  estimation  process  would  have  to  be  reset 
because  every  double  diffeience  ambiguity  would  be 
compromised 

This  process  accouits  for  any  cycle  slip  by  maintaining 
three  independent  RKCPT  solutions  which  uses  the  first, 
secondand  thirdhighest  elevation  SVs  as  the  refo'ence  SV. 
If  the  highest  elevation  SV  loses  lock,  the  final  RKCPT 
solution  simply  implements  another  set  of  ambiguities 
from  another  target  SV.  By  having  three  independent 
RKCPT  solutions  and  ambiguity  estimators,  there  should 
always  be  fixed  ambiguities  available  for  the  final  approach 
solution. 

Independent  Final  Solution.  With  the  high  frequency  of  S  V 
maskings  and  cycle  slips,  a  separation  between  the 
approach  solution  and  the  ambiguity  fixing  process  must 
be  maintained  Tlie  final  solution  only  employs  the 
ambiguities  that  are  known  to  be  fixed,  while  the  Kalmain 
filter  is  resolving  the  ambiguities  that  are  not  fixed 
Therefore  this  pixKess  allows  for  the  highest  amount  of 
ambiguities  to  be  used  in  the  final  solution  always 
anticipating  the  loss  of  SVs  due  to  masking  and  cycle 
slips.  By  allowing  all  ambiguities  that  are  fixed  into  the 
final  solution,  the  availability  of  the  fixed  solution  is 
increased 

The  final  relative  solution  is  simply  a  least  squares  solution 
that  uses  the  sum  of  the  raw  carrier  phase  double  cfiffeiences 
and  either  the  fixed  or  floating  ambiguities  as  range 
measuiements  (see  Equation  3).  This  least  squares  solution 
also  aillows  for  the  needed  high  rate  inatial  data  Inertml 
daita  from  both  ship  and  air  au'e  integcited  in  the  manner 
described  in  the  KCPT  sectiai  to  provide  the  10  Hz 
solution,  and  to  correa  for  GPS  measuiement  processing 
latency.  The  results  in  this  paper  show  the  KCPT  solution 
without  inertial  blending  and  singly  redunchnt. 

RESULTS. 

This  solution  has  been  evaluated  against  two  sets  of  actual 
approach  and  lancfing  data  The  first  set  of  data  was  colleaed 
for  a  Category  in  FAA  contract  in  which  E-Systems 
successliilly  cc^mpleted  100  ^ipproxhes  and  lancfii^s  with  a 
NASA  operated  laser  pnividng  a  truth  referoice.  The 
second  set  of  data  was  colleaed  on  a  helicopter  making 
seveni  jq^pro^dies  to  the  IJSS  ENTERH^ISE  (CVN-65). 

FAA  Category  mb  Data  The  RKCPT  solution  was  used  to 
solve  several  of  the  100  approaches.  Figures  4  and  5  show 
the  position  error  compared  to  an  Ashtech  LI  generated 
solution  known  as  PNAV.  The  Ashtech  KCPT  solution 
was  checked  for  conficbnce  against  a  NASA  operated  laser 
truth  source  The  solid  line  is  the  fixed  solution  error 


truth  source  The  solid  line  is  the  fixed  solution  error 
magnitude  while  the  dashed  line  is  the  OTor  of  the  floating 
solution.  Figures  6  and  7  show  the  ambiguity  error 
compaicd  to  the  true  ambiguity  for  the  best  and  worst 
perfoiming  fixed  base  tpproach  data  We  can  see  that  in 
some  cases  the  Kalman  filter  is  able  to  perfoim  very  well. 
In  the  tuse  of  Figure  6,  all  that  is  requiied  is  to  round  to 
neaiest  ambiguity  to  achieve  the  correct  ambiguity. 
However,  in  the  case  of  Figure  7,  the  ambiguity  neva- 
converges  to  the  correct  solution.  This  is  the  advantage  of 
the  atnbiguity  estimdor.  It  allows  the  more  certain 
ambiguities  to  help  in  the  estimation  of  the  less  certain 
ambiguities. 
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TLSS  ENTER  PR  TSE  Data  This  analyas  used  thirty  minutes 
of  data  colleaed  on  an  SH-60  helicopter  making  four 
e^proaches  to  the  USS  ENTERPRISE  (CVN-65).  Figure  8 
shows  the  ground  trace  of  the  Enterprise  (dashed  line),  and 
the  ground  trace  of  the  heliccpto"  (solid  line).  The  thick 
portion  of  the  helicopter  profile  is  the  time  that  the 
solution  was  not  fixed  For  every  approach,  the  solution 
was  initiriized  to  show  that  the  RKCPT  was  able  to  fix 
ambiguities  in  the  time  allotted  for  a  typied  ^proach. 
Time  to  fix  ambiguities  for  each  approach  was  41,  24,  109 
and  21  second  respectively.  Figures  9  and  10  show  the 
altitude  and  range  plots  with  refeimoe  to  time.  The  dashed 
line  is  high  when  the  ambiguities  are  fixed  Figure  11 
shows  how  the  Kalman  filtCT  and  final  solution  intact. 
The  solid  line  is  the  number  of  SVs  used  in  the  fixed 
solution,  the  dashed  dot  line  is  the  numba:  of  SV 
ambiguities  used  in  the  float  solution,  and  the  dotted  line  is 
the  number  of  SVs  available.  It  shows  how  the  solutbn  is 
constantly  losing  and  resolving  the  ambiguities  to  give  the 
seamless  fixed  ambiguity  opa:alion. 
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CONCLUSION 

This  method  shows  that  is  possible  to  anploy  a  KCPT 
system  aboard  ship  (and  mobile  ground  stations  in  geneid) 
to  piDvitfe  precision  t^proach  and  landing  position 
solutions.  To  operate  an  RKCPT  solution  in  a  shipboard 
environment,  the  process  must  be  able  to  accommodate 
multiple  SV  maskings  and  cycle  slips  on  a  moving 
platfonn.  The  test  data  show  that  the  RKCPT  method  is 
able  to  fix  and  maintain  carrier  phase  ambigiities.  The 
next  stq)  to  take  in  developing  a  GPS -based  shipboard 
landii^  system  is  to  add  the  inertial  blend  filter  described 
above  to  the  RKCPT  solution  to  piovi*  rdative  postion- 
vdodty-acceleaation  data  at  an  adequate  rate  for  established 
precision  approach  efisplay  and  control  algorithms. 
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